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1 Introduction 

The Gamow (or resonance) states are the wave functions of resonances. These 
states are eigensolutions of the Schròdinger equation subject to a "purely outgoing 
boundary condition." The Gamow states were introduced by Gamow [1] in 1928 
to describe a decay (see also [2]). Some years later, in 1939, Siegert made use of 
the Gamow states to obtain a resonance expansion of the scattering function for 
potentials of finite range [3] which was further developed by Humblet and Rosen- 
feld [4]. In 1955, Peierls pointed out the relationship between the residues of the 
propagator at the resonance poles and the Gamow states [5] . The properties of the 
Gamow states and its applications have been considered by many authors, see for 
example [6-20]. A pedestrian introduction to these states can be found in [21]. 

The Gamow states can be used to obtain resonance expansions of wave func- 
tions and propagators. The purpose of this paper is to discuss the main features 
of several of such expansions and to point out some of their differences. For the 
sake of simplicity, we shall restrict our discussion to potentials of finite range. The 
potentials will be three dimensional and sphcrically symmetric except in Section^J 
where they will be assumcd to be one dimensional. 
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2 Berggren's and Berggren-like resonance expansions 

In order to derive Berggren's expansion [8], we start out with a completeness 
relation in terms of the bound and scattering states, 

Nb i-oo 

l = Y^\K n ){K n \+ dk\k+)(+k\, (1) 

n=l J ° 

where 1 is the identity operator, \K n ) are the bound states and \k + ) are the scatter- 
ing states. We are assuming that the potential holds N\> bound states. By deforming 
the integral of Eq. £[J into the contour Tb of Fig-QJ and by using Cauchy's theorem, 
we extract the contributions from the resonance states that are hidden in the con- 
tinuum integral and writc thcm in the same form as the bound state contributions: 



1 = ^ 



K n ){K n \ + Y^\k+){+k n \+ f dk\k+)(+k\, (2) 

n=l n=l •' r B 

where |fc+) is the Gamow state corresponding to the nth resonance, and Nb is the 
number of resonances that lie in between the real axis and the contour Tb- The 
integral in Eq. is called the background. 

Berggren's expansion includes only those few resonances that are supposed to 
be the most important in the energy range of interest. The background integral is 
usually assumed to give a negligible contribution in such energy range. However, the 
background need not always be negligible. For example, virtual poles influence low 
energy scattering, and the background integral in Eq. @ cannot be neglected at low 
energies. (Of course, one can easily deform the contour Tb to include explicitly the 
effect of the virtual poles.) As well, even when the background integral is negligible, 
it is never zero. This is due to the fact that the bound and resonance states do not 
form a complete basis, and one has to include an additional sets of kets to obtain 
a complete basis. When the background integral is not negligible, which is the 
usual case encountered in applications, one needs to estimate it. For Berggren's 
expansion, a discretization of the background integral yields such estimation [16]. 

One can construct Berggren expansions of the Hamiltonian H, the resolvent 
l/(z — H) and the evolution operator e~ lHt by simply letting those operators act 
on Eq. ©. 

Berggren's expansion @, sometimes with slight modifications, is the one most 
oftenly used in nuclear physics [12, 13, 16, 19]. 

Needless to say, the completeness relation (j2J is a formal expression that must 
be understood within the rigged Hilbert space as part of a "sandwich" with well- 
behaved wave functions / and g: 

JVb N B „ 

(f, 9 ) = '£(f\K n )(K n \g) + '£(f\k+)( + k n \g) + / ^k (f\k+) (+%) . (3) 
n=l n=l 

Expression @ is vàlid only when the wave functions f(r) and g(r) fall off at infmity 
faster than any exponcntial. 
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As we mentioned above, Berggren's expansion includes only the resonances that 
carry the most influence in the energy range under consideration. One can incor- 
porate the contribution of othcr resonances by enclosing other poles of the fourth 
quadrant of the fe-plane. For example, by using the contour Tm of Fig. one 
obtains [17] 

N b N M . 

l = J2\Kn){K n \ + J2H){ + k n \+ dfe|fe+)(+fe|, (4) 
n=l n=l "^ r M 

where Nm is the number of (proper) resonances in between the contour Tm and 
the real axis. By substituting |fc + ) = S(k)\k~) into Eq. QJ, where S(k) is the S 
matrix and \k~) is the "out" Lippmann-Schwingcr ket, and by using the contour 
Teg °f Fig. HI onc obtains another expansion [11]: 

l = J2\Kn)(K n \+J2\kn)( + kn\+ dk \k~ ) S (k) {+ k\ , (5) 

n=l n=l 

where N-qg is the number of resonances in the fourth quadrant of the complex 
plane. One can also include virtual states in an obvious way 

Likewise the completeness relation (01, the completeness relations (@J and © 
are to be understood as part of a "sandwich" with well-behaved functions / and g. 
However, unlike the completeness relation @, the completeness relations (0J and 
JSJ do not make sense as they stand. The reason is that, in order to properly derive 
those expansions, it is necessary that the analytic continuation of the integrands 
(g\k + )( + k\f) and (g\k~)S(k)( + k\f) tends to zero as k tends to infinity in the fourth 
quadrant of the complex plane. However, those integrands diverge exponentially in 
that limit. For example, if f(r) is an infinitely differentiable function with compact 
support, then its wave number representation /(fe) = ( + fe|/) diverges exponentially 
on the infinite arc. Therefore, one has to use either a regulator or a time-dependent 
approach: 

N b N BG 

e ~im = Y j e~ iK " t \K n ){K n \+Y J e- ifc " t |fe+)( + fc„|+ / dfe e - ifc2t |fe + )(+fc| , t > 0. 

n=l n=l 

(6) 

The expansions (J2J and © are now understood as the (singular!) limit when t — > 
of the expansion ©. Note that the expansion (JüJ is vàlid for t > only, yet 
another reminder that resonances ought to be understood in a time-asymmetric, 
time-dependent manner. Wc shall further discuss time-dependent expansions in the 
next section. 

A side effect of the divergence of the integrand (g\k~)S(k){ + k\f) at the infinite 
arc is that the so-called "Hardy axiom" [11,22] is flawed. The reason is that the 
"Hardy axiom" assumes that (/|fe~) and ( + fc|<7) are Hardy functions, and therefore 
(/|fe _ ) and ( + k\g) should tend to zero in the infinite arc of the fourth quadrant 
of the complex fc-plane. On the contrary, the quantum mechanical wave functions 
(/|fe _ ) and ( + k\g) diverge (exponentially!) in that limit, and therefore they cannot 
be Hardy functions. 
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3 Time-dependent expansions 
3.1 Expansion involving proper resonance poles 

As is well known, the solution ip(t) to the time-dependent Schròdinger equation 
may be written in terms of the retarded time evolution operator exp(— iHt), where 
t > 0, and of the known arbitrary initial state ip(0) as 

m = e- im m . (?) 

In coordinate representation, the retarded Green function g(r,r';t) — (r\e~ lHt \r') 
may be written, using the Laplace transform, as [9] 

g(r, r'; t) = ± f dfc G+(r, r'; fc) e"^ 4 2fc , (8) 

where G + (r, r' ; fc) denotes the outgoing, time- independent Green function. The 
contour C runs parallel to the positive imaginary fc-axis and bends close to the 
origin to run parallel to the positive real fc-axis, staying always in the first quadrant 
of the complex fc-plane, see Fig.^ It is possible to expand the time-dependent Green 
function in terms of the resonance states plus a background integral by deforming 
appropriately the contour C in the fc-plane [9]. Since the variation of G + (r, r';fc) 
with fc complex is at most exponential [23], the behavior of the integrand with fc 
in Eq. (JSJ is dominated by exp(— ifc 2 í). A convenient choice is to deform G to a 
contour involving two semi-circles C s along the second and fourth quadrants of the 
fc-plane plus a straight line Cl that passes through the origin at 45° off the real 
fc-axis. In doing so, one passes over some poles of the outgoing Green function. In 
general, these include bound states and the subset of complex poles associated with 
the so-called proper resonance states, i.e., those poles for which Refc p > Imfc p . By 
extending the integration contour up to infinity, one obtains that the semi-circles 
C s yield a vanishing contribution, so one is left with a infinite sum of terms plus 
an integral contribution [9] : 

oo 

g(r,r';t) = 5^« p (r)u p (r / )e- I ** í 
P =i 

+ í dfcG + (r,r';fc)e- ifc2t 2fc, (9) 
27r Jc L 

where, without loss of generality and for the sake of simplicity of the expressions, 
we have omitted the bound states. In the last equation, the sum runs through the 
proper resonance poles -hence the subscript p- of the outgoing Green function. In 
deriving expression (jj), one uses that the residue of G + (r, r'; fc) at the complex pole 
fc p is Up(r)u p (r) /2k p , and that the resonance states are normalized according to the 
following condition [9] : 

[ R d rul(r)+i^-=l, (10) 
Jo 
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where R denotes the radial coordinate from which the potential vanishes. Note that 
the integral term in Eq. Q may be written as an integral along the contour Tm- 

3.2 Expansion involving the full set of poles and resonance states 

One may obtain an expansion involving the full set of bound, virtual, resonance and 
anti-resonance states (the latter associated with the third-quadrant poles of the S 
matrix) by noting that the contour in Eq. JHJ) may be deformed into a semi-circle 
C' s along the third quadrant of the fc-plane plus an integral term along the real 
fc-axis. By extending the integration contour up to infinity, we again obtain that 
the contribution of the semi-circle vanishes and we are left with the expression 



i í +0 ° 2 
g(r,r';t) = — dk G + (r, r; fc)e~ lfc 1 2k . (11) 

It turns out that under the condition r, r' < R, one may use the Cauchy expansion 
of the outgoing Green function 

n+l i u\ 1 u n {r)u n (r') 

G + (r,r;fc) = - £ , r,r<R, (12) 

n— — oo 



where the anti-resonance poles correspond to negative integer vàlues of n. The above 
expression still holds whcn eithcr r or r', but not both, equals R. To our knowledge, 
the above result was first reported, for a solvable model, by More [24]. It has been 
proved in the s-wave case for potentials of finite range by García-Calderón and 
Berrondo [25], and it has been used by several authors, see for example [26-28]. The 
expansion given by Eq. 112|) implies that the Gamow states satisfy the relationships 

± !h^n =0 , r y <R , (13) 

fon 

n— — oo 

and 

2 

Substitution of Eq. (O into Eq. JT]]) yields [29] 

OO 

g(r 7 r'-t)= ^ u n {r)u n (r')M(k n ,t) , r,r'<R, (15) 

n— — oo 

where the M function is dehned as 

M(k n ,t) = 1- í 

27T J_, 



^ :v 

- Yl u n (r)u n {r r ) = S(r - r') , r,r' < R. (14) 



n— — oo 



p -i/c 2 í 



= \w{iy n ). (16) 
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The function w is the complex error function [30]: 

w(z) = cxp(— z 2 )erfc(— iz) . (17) 
The argument y n is given by 

Vn - -e-V • (18) 

One may write the above solutions in a form that exhibits explicitly its exponen- 
tial and non-exponential contributions by writing them in terms of an expansion 
involving the complex poles of the fourth quadrant of the fc-plane: 

oo 

g(r, r'; t) = ^{^(r^r 1 )e~ ik l l ~[u p (r)u p (r')M(-k p , t)-u* p {r)u* p {r')M{-k* p , t)}} , 
p=i 

(19) 

where we have used that the poles of the third quadrant satisfy, from time reversal 
considerations, that /c_ p = — k* and that = u*. Also, we have made use of the 
following symmetry of the w functions [30] : 

w(iy p ) = 2é>l - w(-iy p ) , (20) 

which holds when the argument y p lies within the limits 

tt/2 < argy p < 3tt/2. (21) 

This last condition is fulfilled by the proper poles of the fourth quadrant. 

We note that, contrary to all previous expansions, the expansion (|15fl nrst yields 
an estimation of the background integral when r, r 1 < R, second it is vàlid when 
r, r' < R, and third it includes the contribution of the anti-resonance states explic- 
itly. 

4 Resonance expansions and the steepest descent method 

Resonance expansions of the wave function in coordinate or momcntum rep- 
rescntations, and of survival amplitudes arise naturally when calculating the cor- 
responding integrals with the steepest descent technique. The important point is 
that, even though the poles are the same as in other approaches (say poles of the 
transmission amplitude in one-dimensional scattering) , their contributions may dif- 
fer. In this paper, we shall illustrate the bàsic idea in coordinate representation and 
for one-dimensional scattering off a potential with support [0, d], although similar 
manipulations can be performed in other cases. Instead of the wave number k, in 
this section we shall use p/h, p being the momentum. 

Assume that the wave packet is initially confined to the left of the potential so 
that it can be written, in terms of its initial momentum representation 4>(p) and 
the scattering eigenstates of H, tp p , as 

(a#(t)>= / àp(x\%)e- {Et l h <t>(p) (22) 
Jc 
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where 

jpx/h + R(p) e -ipx/h x < o 
T(j) y P x/n x > d ^ 

and C goes from — exo to oo above all singularities. The coeficients R{p) and T(p) 
are the reflection and transmission amplitudes for p > 0, continucd analytically 
elsewhere. The transmitted wave packet may then be written as 

(a#(í)> = h- 1/2 í dpe ixp/h e- iEt/h T(p)<j>(p), (24) 
Jc 

and the contour C may be deformed along the steepest descent path. 
These integrals may be written in the form 



l=í dke-^+^gik). 
Jc 



(25) 



In simple cases, g(k) is a meromorphic function. The saddlc point of the exponent is 
at k = —b/2a, and the steepest descent path is the straight line Im(fc) = — (Re(fc) + 
b/2a). By completing the square, introducing the new variable u, 

u=(k + b/2a)/f, / = (l-i)(m7i/í) 1/2 , (26) 

which is real on the steepest descent path and zero at the saddle point, and mapping 
the contour to the u-plane, the integral takes the form 

' ht f í due- u2 G(u), (27) 



l = e" 

where G(u) = g[k{u)\. It is now useful to separate the pole singularities explicitly 
and write G as 

G W = E^ + ^)' (28) 

3 3 

where Aj/f is the residue of G(u) at u = Uj and the remainder, H(u) 1 is obtained 
by substraction. Note that H(u) is an entire function if G is meromorphic. 

The integral X is thus separated into two integrals, X = 2' + X" . The first one 
may be reduced to known functions by deforming the contour along the steepest 
descent path (real-u axis) and taking proper care of the pole contribution, 



jy = e im X ynty A f du J±_L_ 
j Jc u u-uq 

jm X */ntJ2A jW (- Uj ). (29) 



-i7re" 



The second integral, which involves the remainder H, must in general be evaluated 
numerically, 

I" = e ímx2 l u f ( duH(u). (30) 
Jc„ 
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However, the computational effort is greatly reduced by deforming the contour along 
the steepest descent path, too. When it is an entire function, it may be expressed 
as a series by expanding H(u) around the origin and integrating term by term, 



-r" Amx 2 /M /• 1/2 

m = e jtt ' 



H(u = 0) 



(31) 



E 



1 x 3 x ... x (2n - 1) 
2™(2n) 



H {2n \u = 0) 



In practical applications, the first term alone gives already a very good approx- 
imation. Note the bàsic role of the w-functions, which may be thought of as the 
elementary transient mode propagators of the Schròdinger equation. Long and short 
timc formulae are easily obtained from asymptotic expansions of uu(z). 

For applications of this technique to describe different transient phenomena, see 
for example [31-33]. 



5 Complex, effective Hamiltonians 

As we shall see in this section, truncation of resonance expansions enable us to 
understand the origin of complex, effective Hamiltonians. 

If we apply the Hamiltonian to the Berggren expansion (J2J) , we obtain 

H = Y J K 2 n \K n ){K n \ + Y J kl\k+){+k n \+ dkk 2 \k+){+k\, (32) 

.. — 1 .. — 1 -^Tr 



Thus, the Hamiltonian can be split into a sum of bound, resonance and background 
contributions: 

H -^bound ^resonance ~l~ -^background ■ (^^) 

In matrix notation, one can write this equation as a diagonal matrix, 
/ K 2 \ 



H = 



V 



K 2 



k 2 



k 2 



(fc 2 )bg / 



(34) 



Whcn, for example, the first and second resonances are the most important to the 
problem under consideration, we can neglect everything but the contribution of 
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those two resonances. The corresponding cfïcctive, complex Hamiltonian arises in 
a natural way: 



When double poles come into play, the effective Hamiltonians have non-diagonal 
terms [17]. Note that, contrary to PT-symmetric Hamiltonians, the eigenvalues of 
this effective Hamiltonian are complex. 



RM acknowledges financial support from the Basque Government through rein- 
tegration fellowship No. BCI03.96. 
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